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Abstract 

We investigate a semiclassical conductance for ballistic open three- 
dimensional (3-d) billiards. For partially or completely broken-ergodic 3-d bil- 
liards such as SO (2) symmetric billiards, the dependence of the conductance 
on the Fermi wavenumber is dramatically changed by the lead orientation. 
Application of a symmetry-breaking weak magnetic field brings about mixed 
phase-space structures of 3-d billiards which ensures a novel Arnold diffusion 
that cannot be seen in 2-d billiards. In contrast to the 2-d case, the anomalous 
increment of the conductance should inevitably include a contribution arising 
from Arnold diffusion as well as a weak localization correction. Discussions 

are devoted to the physical condition for observing this phenomenon. 
PACS numbers: 03.65Sq; 05.45.+b; 73.20.Dx. 
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I. INTRODUCTION 



Concave and convex billiards as prototypes of conservative chaotic systems have received 
a growing theoretical and experimental interest in the fields of nonlinear dynamics and sta- 
tistical mechanics [1]. These billiards are nowadays fabricated as quantum dots or antidots 
in ballistic microstructures where the system size is much less than the mean free path 
20/xm) and larger than the Fermi wavelength (A ~ 50nm)[2]. Since the striking exper- 
iments on the magneto-conductance of nanoscale stadium and circle billiards with a pair 
of conducting leads [2], quantum transport phenomena have been intensively investigated 
[3]. Some salient features of experiments well agreed with the semiclassical theory, where 
a semiclassical conductance formula was derived for two-dimensional (2-d) chaotic billiards 
by combining the semiclassical Green function with the Landauer formula [4,5]. Subse- 
quently, the semiclassical conductance formula was extended for integrable 2-d biUiards [6] , 
and detailed numerical computation on circle billiards indicated that the conductance also 
depends on the lead orientation [7]. Recently, the dependence of the conductance on the 
lead orientation was found experimentally in a ballistic square billiard [8]. 

At present, another important object of research on quantum transport is the three- 
dimensional (3-d) biUiard [9]. The ballistic 3-d quantum dots will be fabricated, for in- 
stance, by exploiting drying etching processess with focused-ion beams applied to AlxGai_x 
As/GaAs/AlxGai_xAs double hetero-structures, which will be less than both the elastic 
mean free path and phase coherence length. The subject of this paper is to extend the 
semiclassical conductance formula for 2-d billiards to 3-d billiards and to investigate the 
quantum transport in open 3-d billiards from a viewpoint of the Arnold diffusion, a key 
concept in high-dimensional nonlinear dynamical systems. While Arnold diffusion is a clas- 
sical phenomenon, an anomalous phenomenon like this in classical mechanics should have 
a quantum counterpart, which can only be analytically revealed by using the semiclassical 
theory. 
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II. SEMICLASSICAL CONDUCTANCE FOR OPEN 3-DIMENSIONAL 

BILLIARDS 



The conductance for open 3-d billiards connected to a pair of rectangular parallel- 
piped lead wires can be expressed by the transmission coefficient T as G — — 
^ Sn,m=i l^nml^, where tnm IS the S-matrix element between the incoming mode m = 
{mx,my) and the outgoing mode n = {nx,ny), and the double summation is taken over 
all propagating modes. The transmission amplitude connecting those modes at the Fermi 
energy Ep reads as [10,11] 

j dx' dy' dx dy %l)*n{x'\y'')G{x\y'\ z \x ,y , z\ EF)'il^m{x\ y'), (1) 

where and are the longitudinal velocities of electrons for the incoming and outgoing 
modes, respectively, x' ,y' , z and x" , y" , z" are the local coordinates for the transverse {x, y) 
and longitudinal [z) directions (inward the billiard for z and outward for z") of the incoming 
and outgoing leads, respectively, and G{x" ,y" , z" ]x' ,y' , z ; Ep) is the Green function for 
an electron propagating from the entrance to the exit, i^mix'.y) and ipnix'^y") are the 
transverse components of the wave functions at the leads. To simplify the problem, the 
cross section of the leads is assumed to be a square with side length I. Then '^^{x ,y') — 
jsin(mx7cx' /I) sm(my7ry' /I) and ipn{x" , y" ) = y sin(na,7ra;"//) sin(n^7ri/"//). In the ballistic 
regime, the Green function from entrance (r ) to exit (r") can be replaced by the Gutzwiller 
semiclassical path-integral expression [12] 

G{x ,y ,z ;x ,y ,z;Ef)^ (27ri^)2^ ^ | det(-^)|^/' exp(-5,(r , r , - i7ri/,/2), 

^ ' s(r',r") 

(2) 

where Ss and Ug are the action and Maslov index for a trajectory s, and vp is the Fermi 
velocity. The above semiclassical Green function was widely apphed to derive the level 
density of both chaotic and regular billiards [12,13]. Here, to avoid direct transmission 
without any bouncing at the billiard wall we may place a stopper inside the billiard. 
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In the integration in Eq.(l) we introduce the entrance angle Q' between the initial di- 
rection of the orbit and the z axis, and the angle between the projection of the initial 
direction onto the x —y plane and the x axis. Angles 9" and Lp" are defined for the exit in the 
same way. Then, we can rewrite | det(^p7f )| for each trajectory s as | det(^^7^^;^)| = p^p-^s 
where Ag = I '^'"f// — l^-f^)!- For large mode numbers, the integration over the 

' cos ^ ax oy oy ax ' ^ ^ ' ° 

transverse coordinates in Eq.(l) can be performed in the stationary-phase approximation, 
and the transmission amplitude becomes 

tnm^-ifi^ E sgn(m^)sgn(mj,)sgn(nJsgn(nj,)(cos^'cos^'V/2(^/5,)^/2gt5.-iPW2 

■P^ s{mx,fhy,nx,ny) 

(3) 

with 

^ ir . ^ II , do" do" II dip" dip" ^11. „ ^11 , dp" do" , n dO" dp" , 
= sm2(^ cos'e -^^siv?e )+sm20 {^—cos^cp - „ /„ sm' cp 

oy ox ay ox oy Ox oy Ox 

r. ^ I , do' do' ^1 dp dp . n I. . ^„i ,dp' do' o / dO' dp . n /,,, 

where fhx — ±mx,fhy — i^ruy, fix — ±nx,ny — ±%, and Sg — S{r" r'^, Ep) + hnfhxx'g/l + 
hnifiyy'g/l — hnfixx'g/l — hnfiyy'g/l, and z> is the Maslov index which includes an extra phase 
coming from the possible sign change in each of the four Fresnel integrals. The above 
result means that only those isolated trajectories connecting a pair of transverse planes at 
discrete angles sin O' cos (p' — ffix'K /k pi, sin O' sin p) — fhyn /kpl, sin O" cos p^" — nxn/kpl, and 
sin O" sin p" = Uyii / kpl dominate the conductance. 

Using Eqs.(3,4), we shall first analyze a completely chaotic 3-d billiard with a pair of 
conducting leads. In this case, the electron injected from the incoming lead will bounce at 
the billiard boundary in an ergodic way before reaching the exit, leading to amplitudes of the 
same order for both transmission and reflection coefficients. In the large mode number case, 
the summations over modes are replaced by integrations and the transmission coefficient 
T — J2n,m l^nmP IS rewritten as 

T = n^'^i. V Vcos^'cos^"(^)^/2et(^-^")-f(^-'^") 

fF mx,my,nx,ny s,u ^s^u 
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/ sm6''(i(sin6'') / li^^' / sin^"(i(sm6'") / (iy?" cos O' cos ^" 

v/ — 1 J— TT J— 1 J— TT o „ 



Noting the dimension of (^^)^/2 to be t^, and using the diagonal approximation, the above 
integration can be evaluated as T oc {kpiy. The conductance for 3-d chaotic billiards is 
thus given by G oc ^{^piy without any dependence on the lead orientations. 

III. COMPLETELY OR PARTIALLY BROKEN-ERGODIC 3-D BILLIARDS 

However, the situation will be dramatically changed for the case of completely or par- 
tially broken-ergodic 3-d biUiards. As examples, we choose S0(2)-symmetric biUiards: A 
completely integrable (: broken-ergodic) 3-d billiard is available by rotating, e.g., the 2-d 
ellipse billiard (on the Y — Z plane) around the Z axis (see Fig.l). Note: X,Y,Z are the 
global coordinates; A partially broken-ergodic 3-d billiard is also available in a similar way 
by rotating, e.g., the 2-d stadium. All these S0(2)-symmetric 3-d billiards have the an- 
gular momentum Lz as a constant of motion. In an open system version of each of these 
billiards, one may place in its inside a smaller stopper of the similar shape, to avoid direct 
transmission of electrons. The resultant 3-d shell biUiards still retain the S0(2) symmetry. 

Firstly, consider the vertical case (i) when the incoming and outgoing leads are connected 
with the billiard vertically at points A and C, respectively (see Fig.l). Electrons to reach the 
exit C should have the vanishing angular momentum, Lz — 0. Therefore, only the electrons 
with the initial velocity vector at the entrance A lying in the Y — Z plane can reach C 
because of the conservation of Lz; Other incoming electrons falling into the trajectories out 
of this plane cannot reach C and should return to A. In deriving the transmission coefficient 
for this lead orientation, the integration prior to Eq.(5) is reduced to the line integration 
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FIGURES 




FIG. 1. Schematic ihustration of the 3-d elliptic shell billiard. Dual 2-d ellipses sharing the 
same foci inY — Z plane are rotated around the Z axis to generate a 3-d elliptic shell billiard. The 
longer radii of outer and inner ellipses are taken as 23f? and 1.855R. The shorter radius of the outer 
ellipse is 0.85R. The side length / of the square lead is taken as 5R/20. 

over 9' and O" only, yielding the conductance G oc ^kpl (linear in kp^-)- (The reflection 

coefficient, which is obtained by applying the same procedure as in the chaotic billiard, is 
proportional to {kplY .) We shall then consider the parallel case (ii) when the incoming 
and outgoing leads are parallelly placed at points Q and C, respectively (see Fig.l). Each 
electron incoming from the entrance Q has the vanishing angular momentum, Lz = 0, and 
thereby can reach C or return to Q with probability of the same order. Therefore, for this 
lead orientation, G oc ^{kpl)'^- 

Generally, in completely or partially broken-ergodic 3-d billiards, we have G oc ^{kpl)^ 
with < 7 < 2, and the exponent 7 is determined by how the billiard symmetry is broken 

6 



by the lead orientation. It should be noted: open 2-d billiards always share G oc ^kpl [4], 
independently of both integrability and lead orientations. 

IV. EFFECTS OF SYMMETRY-BREAKING WEAK MAGNETIC FIELD 

The magneto-resistance is also a very important problem in 3-d billiards. It is well 
known that in ballistic 2-d billiards the decrease of resistance with increasing magnetic field 
(: negative magneto-resistance) occurs due to the weak localization based on the quantum 
interference between a pair of time-reversal symmetric orbits [2-5]. What kind of additional 
novel phenomena will be expected in S0(2)-symmetric 3-d billiards with the vertical lead 
{A, C) orientation when a symmetry-breaking weak magnetic field B will be applied along 
the X axis? In contrast to phase-space structures in 2-d systems where each of chaotic 
zones generated by weak B field is mutually separated by the Kolmogorov-Arnold-Moser 
(KAM) tori, those in 3-d systems consist of chaotic zones mutually connected via narrow 
chaotic channels called the Arnold web (AW) [14-16]. In this case, electrons incoming in an 
arbitrary direction, which permeate (chaotic) zone by zone in an ergodic way through the 
AW, can show a global diffusion, i.e., Arnold diffusion (AD). Therefore, a B field increases 
the transmission channels and the number of stationary points in the integration in Eq.(3), 
and the resultant integration in Eq.(5) over tp and tp" as well as d' and d" indicates G oc 
^{kpl)^- In general we shall expect a more general result, G oc {kplY with 1 < /3 < 2, 
and the exponent (3 depends on the width of AW. A crucial point is that a B field has 
increased the exponent in the power-law behavior of the semiclassical conductance in case 
of the vertical lead orientation. This anomalous phenomenon is just the negative magneto- 
resistance beyond the weak localization. We here keep to employ the terminology of the 
negative magneto — resistance wherever a B field reduces the zero-field resistance. 

On the other hand, for the parallel lead (Q, C) orientation where each electronic trajec- 
tory lies almost in the vertical plane, G oc ^{kpiy is retained, irrespective of the absence or 
presence of B field; We can expect only the B field-induced increment of the proportionality 
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constant, which imphes a conventional negative magneto- resistance caused by suppression 
of the weak locahzation. 

Before moving to the next section, we should note: Originally, the negative magneto- 
resistance was evidenced in the numerical analysis of the conductance of 2-d quantum dots 
[4]. The accompanying theory, based on the resistance, manipulated backscattering orbits 
and attributed this phenomenon to the weak localization [4]. Later works [5] based on the 
conductance, which introduced the small-angle-induced diffraction effect together with an 
idea of interference between a pair of partially time-reversal symmetric orbits [5], swept 
away a cloud (i.e., breakdown of unitarity) hanging over the above theory. All our results 
in the present section have been obtained just within the same primitive analysis of the 
conductance as the original one for 2-d dots. 

V. NUMERICAL RESULTS FOR OPEN 3-D ELLIPTIC BILLIARD 

To numerically verify the above prediction, we shall choose a S0(2)-symmetric 3-d el- 
liptic shell billiard with the vertical lead orientation. This ellipsoid is a typical example 
of completely integrable S0(2)-symmetric 3-d billards. The system in the presence of a 
symmetry-breaking field B parallel to X axis will be examined. To understand the classical 
dynamics of electrons, we analyze as a Poincare surface of section the X — Z plane with 
longitude S and velocity component vz/v^ as Birkhoff coordinates [17]. {v^ = is the ve- 
locity unit corresponding to the lowest mode.) For B — 0, the trajectory from the entrance 
A with an arbitrary initial velocity vector lying out oi Z — Y plane proves to be confined to 
a torus (see Fig. 2(a)), failing to reach the exit C, i.e., S = tt/2. Let r] (= ^/^c) be the ratio 
between the cyclotron radius 3?c and the characteristic length 3? of the cross-sectional ellipse 
(see the caption of Fig.l). Then B field can be expressed as B — r]Bo with Bq — rngV^/e^ 
{Bq corresponds to 3?). In case of a weak field {B = 0.02Bo), Poincare sections for the 
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FIG. 2. Poincare sections {X — Z plane) for the trajectories emanating from point A. Lon- 
gitude S (— 7r/2 < 5 < vr/2) and velocity component vzlv^ [y^ = j^) are chosen as Birkhoff 
coordinates. S = and S = ir/2 imply entrance A and exit C, respectively, (a) B = and ini- 
tial velocity vx/v'^ = 0.821, fy/fO = 0.549, = 0.157; (b),(c) B = OMBq and the same 
initial velocity as (a): (b)P = 500; (c)P = 17000. (d)-(f) B = 0.02So and initial velocity 
vx/v° = 0M3,vy/v^ = 0.538, vz/v° = 0: (d)P = 500; (e)P = 4100; (f)P = 16000. P is the 
number of bouncings at the outer billiard intersecting X — Z plane. For the meaning of Bq, see 
the text. 

trajectory with the same initial velocity vector as in Fig.2(a) are given in Figs. 2(b) and (c) 
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for increasing numbers of bouncings at the outer billiard wall, and similarly those for the 
trajectory with the initial velocity vector lying in the X — Y plane are given in Figs.2(d)-(f). 
One can find: as the B field is switched on, the orbit trajectory first leaves the initial torus, 
enters an outer layer, and, repeating a similar process, diffuses over more and more distant 
layers, exhibiting a phenomenon of the Arnold diffusion (AD). Even though a symmetry- 
breaking field is weak enough to keep the orbit almost straight, AD makes it possible for 
the trajectories with initial angular momenta Lz 7^ to reach the exit C, i.e., S = 7c/2 (see 
Figs. 2(c) and (f)) and eventually to contribute to conductance. 




Kpl/Tl 



FIG. 3. Numerical results for semiclassical conductance versus Fermi wavenumber for A,C 
leads orientation. Circle and square symbols are for B = 0.02i?o and B = 0, respectively. Solid 
quadratic and dashed linear lines are the corresponding fitted curves for their means. 
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To proceed to the numerical calculation of semiclassical transmission amplitudes in 
Eqs.(3,4) for the 3-d eUiptic biUiard with the vertical lead orientation, let us define the 



4x4 monodromy matrix M as 



= M 



or 



. At the entrance we introduce a pair 



of mutually orthogonal coordinates r]_ and rj_ lying in the transverse plane perpendicular 
to the initial direction of each orbit. Similarly r]_ and r\ are chosen at the exit. Using 
elements (Mjj) of M matrix, we rewrite the necessary factors for the numerical computation. 
For example. ^ = "S;^:;;!;' = ^[M„(co3/cos9"M,, - sin^M,) - 

M42(cos/cos^"Mn - sin v9"M2i)]/(MnM22 - M12M21) and |^ = .^i'/^k ^ = 

( sec ip" CSC 6" /Mi4 + esc ip" /M24). While M and u (Maslov index) can be computed 
by the method suggested in Refs.[9,18], the first and last rotation angles in M should be 
determined by using r\ , and r]_ , rj_ chosen as above, since, after folding the orbit, 
the initial and final coordinates should become identical. Combining all these factors with 
Eqs.(3,4), the semiclassical conductance is calculated explicitly. The result, which includes 
both the classical contribution and the quantum correction in an inseparable way, is plotted 
in Fig. 3. Although the data show small fluctuations, the fitted curves for their means are 
described by G oc ^{^fI)"^ and G oc ^kpl in cases of B — 0.02So and B — 0, respectively, 
confirming the analytic issue suggested in the previous section. 

The increase of the exponent in the kp dependence of G upon switching on a weak field, 
which cannot be seen in 2-d cases, provides a numerical evidence that the negative magneto- 
resistance (: a difference between two hues in Fig. 3) in 3-d systems comes from AD as well 
as the weak-localization correction. 



VI. SUMMARY AND DISCUSSIONS 

We have investigated semiclassical quantum transport in three-dimensional (3-d) baUistic 
quantum billiards, and explored novel phenomena that cannot be seen in 2-d billiards. For 
partially or completely broken-ergodic 3-d biUiards such as SO (2) symmetric biUiards, the 
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dependence of the conductance G on the Fermi wavenumber /c^? is dramatically changed by 
the lead orientation. In case of the vertical orientation, a symmetry-breaking weak magnetic 
field turns out to increase the exponent vnkp dependence of G. In marked contrast to the 2- 
d case, the anomalous increase of the conductance should be caused by the Arnold diffusion 
(AD) as well as a weak locahzation correction. 

We should note the following: (i) In an extremely weak field case when the width of Arnold 
web (AW) is thin enough, we shall see G oc {kplY with 1 < /3 < 2. Our numerical calculation 
in the previous section suggests that, to have /3 = 2 in a very large kp regime, 3?/3?c should 
be no less than 10~^; (ii) role of a stopper is also essential. Reflection at an inner convex 
wall should increase the instability of orbits, which is favourable for a genesis of AD; (iii) 
to clearly observe AD, the symmetry- breaking field should destroy all constants of motion 
except for the energy. From this viewpoint, it is more advantageous to choose a S0(2)- 
symmetric 3-d chaotic shell biUiard generated from a fully-chaotic 2-d biUiard with a hollow 
by its rotation around the Z axis. The billiards of this kind, which are partially broken- 
ergodic, are experimentally more accessible than a (completely broken-ergodic) elliptic shell 
biUiard, and all conclusions for the latter should hold for the former as well. 

Finally, we point out the role of temperature in observing the effect of AD on the con- 
ductance. At finite temperatures, the conductance is related to the scattering matrix as 
G — ^ /f^ (ie[— Y.a,in l^n,m|^ with / the Fermi distribution function. Then the integra- 
tion over energy will give a temperature factor smh{nkj^L/lhvF)) '^^^^ ^ the length of trajecto- 
ries. For the vertical lead orientation, this factor suppresses very long trajectories reaching 
the outgoing lead along the Arnold web, and thereby the negative magneto-resistance is 
described only by the weak-localization correction. On the other hand, if the temperature is 
lowered, the self-averaging over kp is not permissible, leading to a suppression of the weak- 
localization correction. Then the increment of the magneto-resistance is mainly caused by 
the effect corresponding to the Arnold diffusion. We hope that a future experiment on 3-d 
ballistic cavities can verify our prediction. 
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